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Scene Matching on Imagery 


e There are a plethora of algorithms in existence for automatic 
scene matching, each with particular strengths and weaknesses 


e SAR scenic matching for interferometry applications is difficult 
for several reasons: 


— Thermal noise causes the fine structure of two SAR images to 
differ, adding noise to the correlation measurements. 


— Geometric speckle noise is similar in images acquired in nearly the 
same geometry, but as the interferometric baseline increases, 
differing speckle noise corrupts the matching correlation. 


— Scene decorrelation is another form of speckle noise difference 
that corrupts matching correlation. 


e SAR scenic matching for mosaicking applications involves 
greater challenges, including severe speckle noise differences, 
layover and shadow effects. 
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Automatic Scene Matching 


Find overlap region and sample points at specified spacing in along 
track and cross track direction. 


Typical window sizes are 64x64 pixels for image data and 64x64 or 
128x128 pixels for height data. 


Uses a modified Frankot's method for rejecting bad matches and to 
provide an estimate of the match covariance matrix. 


Cross correlation uses a normalized mean correlation function, mean 
of search window is calculated over the region in common with the 
reference window. 


Mean of reference 
is constant 


Mean of search 
varies with position 
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Match Correlation Estimate 


Correlation is computed in the spatial domain using a normalized cross 
correlation algorithm. 


Let /,(x) be the image values at a point x in the reference window in 
the first data set, L , the mean of the intensities in that window. 


Let 3 (¥ +d) be the image value at pointž+d in the search window 
of the second data set, and 42 (d) the mean of the intensities in that 
window. 


Viewing each image as a vector in an n-dimensional vector space then 
the cross correlation is computed as 
Fi (h-hb@M-h@) h-hbh@-h@) 
z -7 Cd) = D (d) OT 1G 


where 73 and 71d) are the standard deviations of the image 
intensities in data sets 1 and 2, respectively. 


Space-domain correlation used for speed and because large 
irregular data gaps are trouble free. 
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Baseline Estimation from Offsets 
Baseline estimation in repeat pass interferometry is a rich subject, and 
will be covered completely in a separate module. Briefly: 


Offset field determined from scene matching carries information 
sufficient to reconstruct the baseline to an accuracy com-mensurate 
with ability to co-register images. 


For nearly parallel, smoothly varying, orbit tracks, the baseline can be 
modeled as a simple function of the range P 


B (p) ~ Bsin(9,(p) — a) = App) 
B (p) = Bcos(0,(p) - a) = EE (p) 


è Solve two equations for two unknowns: B, dl 0 


e Evaluating baseline at several widely-spaced along-track 
(azimuth) locations gives azimuth history of the baseline. 
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Converging Baselines 


For orbit tracks with substantial convergence, azimuth offsets 
strongly characterize convergence rates. 
Observation vector: 


{: sag Ap(p;,a, ), Aal p;,a; ), AP( P.A) Aa( P; 14; 41)> En +} 


OB. OB, 
Parameter vector: {B > B.o; Bios We —} 
ds OS 


where y is an image scaling factor resulting from velocity differences 
between the tracks, and (s,c,h) are local coordinates defined by the 
orbit track (see notes on coordinate systems). 
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Baseline Accuracy from Offsets 


e Baselines determined from offset fields are typically accurate to 
a fraction of the scene matching accuracy that depends on the 
model function and the number of offset estimates. Example: 


— Zero baseline image pair with 10 match points across range. 


— Scene matching accuracy of 1/20 pixel resolution at typical single 
look resolution of 15 meters: 75 cm / 3 = 25 cm 


e This accuracy is insufficient for topographic mapping 
applications 


O, = sin cos QO g, 
=320m ; p=830 km, B= 250 m, 0 = 23°, a =0, 0; = 25 cm 


These geometric parameters are typical of ERS with baselines 
suitable for topography. Generally, ground control points tied to 
the unwrapped interferometric phase are required for mapping. 
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Interferogram Formation 





A interferogram 
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Registration Implementation 


e In resampling single look complex image to register properly 
with the reference image, care must be taken in interpolation of 
complex data 

— The azimuth spectrum of squinted SAR data is centered at the 
Doppler centroid frequency - a band-pass signal 

— Simple interpolators, such as linear or quadratic interpolators, are 
low-pass filters and can destroy band-pass data characteristics 

e Band-pass interpolators or spectral methods preserve phase 
fidelity 
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Repeat Pass Interferometry 
Back-end Processing 


Performance Summary: 
L-band: good over most terrain 


C-band: poor over most terrain 
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N-pass Differential Processing 


Short Baseline 
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